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Abstract. 10:1 is a new instrument that has recently been commissioned for the Liverpool Telescope, extending 
current imaging capabilities beyond the optical and into the near infrared. Cost has been minimised by use of a previ¬ 
ously decommissioned instrument’s cryostat as the base for a prototype and retrofitting it with Teledyne’s 1.7pm cutoff 
Hawaii-2RG HgCdTe detector, SIDECAR ASIC controller and JADE2 interface card. In this paper, the mechanical, 
electronic and cryogenic aspects of the cryostat retrofitting process will be reviewed together with a description of the 
software/hardware setup. This is followed by a discussion of the results derived from characterisation tests, including 
measurements of read noise, conversion gain, full well depth and linearity. The paper closes with a brief overview of 
the autonomous data reduction process and the presentation of results from photometric testing conducted on on-sky, 
pipeline processed data. 
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1 Introduction 

The Liverpool TelesoopeP0 (LT) is a fully robotie 2m teleseope sited on the Canary Island of La 
Palma. It is an unmanned faeility insomuch that it does not require an observer during operations. 
Instead, observations are queued by an autonomous scheduler. The scheduler provides the tele¬ 
scope with a ranked list of suitable observing programmes based on current and predicted future 
meteorological data, programme priority and other queue optimisation conditions. A scheduled 
observation can be overriden at any time by a Target of Opportunity (ToO), making the LT partic¬ 
ularly well suited to time-domain observational astronomy where its ability for rapid response has 
been a key asset. 

The LT is a multi-purpose facility hosting up to nine instrument simultaneously. The current 
instrument suite includes a lO'xlO' optical imager (10:0), an integral-field spectrographPI^ (FRO- 
DOSpec), a fast-readout cameraP (RISE), a polarimetei® (RING03) and a low resolution, high 
throughput spectrographP (SPRAT), all of which are available for use on any single night. 10:1 
extends the LT’s imaging capabilities up to a near infrared (NIR) wavelength of ~ 1.7/im where the 
instrument will enable science to be delivered in, although not limited to, the following key areas: 

• Supernovae type la (SNe la). Measurement of NIR light curves for nearby SNe la will 
help to improve our understanding of the physics of SNe la explosions, as well as making a 
contribution to increasing the accuracy of SNe la as cosmological distance measures. 

• Gamma-Ray Bursts (GRBs). Early time followup of the highest redshift GRBs (z > 6), 
including measurements of line of sight dust extinction, will allow for better determination 
of the physical parameters of these events. 

*http://telescope.livjm.ac.uk/ 
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• Active Galactic Nuclei (AGN). Investigations into the variability of the various elasses of 
AGN will allow us to probe the central engines of these objects, helping to discriminate 
between competing models. 

• Other areas of both time-domain (e.g. pixel microlensing, variability surveys) and non time- 
domain (e.g. metallicity gradients in galaxies and followup up X-ray cluster surveys) astron¬ 
omy, proposals for which have already been solicited and undertaken. 

2 Instrument Overview 

10:1 is positioned at the Cassegrain focus of the telescope where, without additional powered 
optics, the plate scale delivered at the focal plane is 97/im/" (2m diameter mirror, f/10 beam focal 
ratio). With a pixel pitch of 18/im, the total FOV of the 2048x2048 pixel detector is 6.27'x6.27'. 
The quantum efficiency of the detector is over 50% between 0.8 and 1.72 /im (peak 86% at 1.5/im), 
although the usable range within this window is limited by atmospheric transmission (see Figure 

ffl). 

The detector itself is a HgCdTe Hawaii-2RG (H2RG) and is controlled via the SIDECAR ASIC 
controller and JADE2 interface card, all provided by Teledyne. Connected to the cryostat are two 
closed-cycle cooling units (ICC Polycold Cryotigers), each charged with PT13 gas, providing a 
total cooling capacity of lOW at a minimum temperature of 80K. Inside the cryostat are two cor¬ 
responding cold heads. One of these cold heads is connected to the detector/SIDECAR mounting 
block. The other is attached to the floor of the radiation shield housing, shifting the bulk of the 
300K blackbody radiative load incident from the warm cryostat chassis onto this thermal path. 
By thermally isolating the two paths, the majority of the cooling capacity provided by a single 
Cryotiger is available to cool the detector and SIDECAR only. 

Both J and H single filters have been procured, although the instrument has no filter wheel to 
allow robotic selection between the two. Eollowing the requirements outlined by science proposals 
already received for this instrument, the H band filter has been preferred over the J. Changing be¬ 
tween these various configurations is not a user selectable option and would require the instrument 
to be taken off telescope. As a possible future alternative, 3'x6' J and H filters have also been 
procured, giving the option to split the field-of-view (EOV) of the detector across the filters. With 
an appropriate dithering strategy, the split configuration would allow the user to observe with both 
filters in a single observation sequence at the cost of reduced field size. Eor example, concurrent 
J and H observations would be obtainable by alternately offsetting the telescope so that when the 
target is being observed in one band, sky is simultaneously being observed in the other. 

The wavelength coverage for each of the filters is shown in Eigure [B Elux in the H band 
will be partially attenuated by the intrinsic 1.72/im cutoff wavelength of the detector. This cutoff 
was chosen to restrict the additional complexity required when observing further into the infrared, 
especially with warm optics and non-IR optimised baffling. 

2 .1 Mechanical 

Eigure[2]shows the prototype cryostat, recycled from the ET’s previous infrared instrument (SupIR- 
Cam) with only minor modifications. The H2RG detector and SIDECAR are housed inside an alu¬ 
minium box which acts as a radiation shield. They are both are mounted upon an invar block, used 
for its low coefficient of thermal expansion, which in turn is screwed to the base of the box by four 
thermally insulating Tufnol supports (see Eigure [3]). This setup decouples the thermal paths from 
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Fig 1 Single and split filter transmittance with modelled atmospheric window derived from ATRANP Top: J band. 
Bottom: H band. The detector cutoff, defined as the point at which the QE drops below 50%, is shown for the H band. 

the two Cryotiger cold heads. The lid of the box holds a 58.5x58.5x5.5mm filter/blank central to 
the optical axis. 

This design closely follows that of RATIR,l^in that instead of potting the 15" flex cable connect¬ 
ing the (cold) SIDECAR and the (warm) JADE2 through the cryostat wall, the JADE2 is instead 
kept under vacuum inside the cryostat. In this configuration, both the power and data lines are 
connected to to a hermetically sealed MIE-Spec socket, with the data transferred externally over 
a multimode OMl fibre optic cable to the instrument control computers (ICCs) located off the 
telescope. 
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Fig 2 Left: Dewar with bottom and top cover removed showing the detector and SIDECAR. The enclosure housing 
the JADE2 is located towards the top of the figure. Right: Mechanical drawing of the cryostat with lid/floor and filter 
box lid removed. The JADE2 is mounted upside-down so that the extent of the resulting extrusion from the main 
cryostat chassis is reduced. 



Fig 3 Left: Detector/SIDECAR and mounting block. The temperature sensor is visible towards the top left, and 
is wound around a spool acting as a thermal anchor. Right: Schematic drawing of the mounting block with features 
labelled. The adjustable mounting bracket allows the SIDECAR to be moved towards/away from the detector, allowing 
some degree of flexibility when fitting the 2" flex cable between the two. The section of the block upon which the 
temperature sensor and resistor are mounted is connected via a copper block to one of the Cryotiger cold heads. 

2.1.1 Additional Bajfling 

Immediately obvious in initial testing was a gross light leak, visible as a >1 kADU/s gradient 
across the array in up-the-ramp (UTR) sequences taken with a cooled blank installed. In UTR 
mode, the array is read out non-destructively at successive integration times throughout an expo¬ 
sure. Upon further examination, this gradient was thought to arise from i) light incident obliquely 
on the cryostat window, and/or ii) thermal/LED emission from the JADE2 card. To resolve the is¬ 
sue, a baffle tube was added to the underside of the radiation shield lid, as well as adding two plates 
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to lightly clamp the 15" flex cable conneeting the SIDECAR to the JADE2 as it feeds through the 
radiation shield. The baffle tube and corresponding light path to the detector are shown diagram- 
matically in EigurelH The telescope uses only optical baffles, and has not been optimised for the 
NIR. 


Upper Baffle 



Fig 4 Section view of the baffle. Light enters from the top. The detector only “sees” the warm science fold mirror, 
located at a distance of 570mm from the detector surface. Units, where unspecified, are in mm. 

To quantitatively assess if the problem had been resolved, the rate of eharge aeeumulation was 
plotted for eaeh pixel in an UTR sequenee following these modifieations. A median rate of <1 
ADU/s was ealeulated following the applieation of a four parameter fit diseussed in ^3.1.11 The 
full frame result is shown in Eigure[5l Aside from the square imprint of the filter housing around 
the edge of the frame, no signifieant variation in gradient aeross the ehip is observed. 

2.2 Cryogenic Considerations and Temperature Control 

In order to maintain temperature, SupIRCam had to be placed on the vaeuum pump at intervals of 
between 2-4 weeks. This was thought to be due to the positioning of the getter ehassis on the floor 
of the radiation shield where the temperature was found to plateau at around 160K. In order to 
make the pumping aetion of the getter more effeetive, it was required to attain a lower temperature 
as the residenee time for a moleeule adsorbed on the surface of a material, r, is approximately 
given byP 


1 

r = —exp 



( 1 ) 


where z/q is the frequeney of oscillation of a molecule on the surface (taken as 10 ^^s) and is 
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Fig 5 2D map of the rate of charge accumulation for each pixel in an UTR sequence. The panels on the right show 
the median rate along both collapsed axes. The effect of different DC offsets for each output channel is evident in the 
topmost right panel. Units are in ADU/s. 


the activation energy of desorption. From Ref. [TOl an ^ 40kJ/mol is assumed (typical of 
physisorption binding energies), giving a residence time of the order of seconds at 160K. After 
repositioning the getter directly on the detector cold head, the temperature of the getter material 
now reaches about 90K. Consequently, the residence time of adsorbed molecules, and thus the vac¬ 
uum hold time, should have increased significantly. In practice, however, other factors such as the 
initial conditions of the getter (and consequently adsorption capacity) make this time an unachiev¬ 
able upper limit of the actual vacuum hold time. The cryostat has been empirically found to have 
a hold time of about 8 months, although recent testing of a strategy to purge the getter material 
through baking and saturation with nitrogen is yielding promising results for future maintenance. 

In order to servo the detector temperature, a LakeShore DT-47I-CO silicon diode has been 
clamped to the detector mount (see Figure [3]) alongside a single 25^1 resistor used as a heating 
element. Both the diode and resistor are connected to a Lakeshore model 331 controller, with the 
diode having both load and sensing lines. This configuration provides a stability of 4.2-5.5mK 
over a temperature range of 77-300K, with the resistor capable of outputting between 0.25W and 
25W of heater power when connected to the LakeShore’s primary control loop. To mitigate the 
possibility of exceeding a rate of temperature change greater than the recommended IK/min, the 
heater power has been set to a maximum of 2.5W. A separate DT-471 built in to the SIDECAR 
is also connected to this controller to monitor the ASIC temperature, but is not connected to a 
control loop. As this generation of controller has only a serial interface, a custom transceiver with 
a built-in ethernet server is used to control it over the telescope network. 

2.3 Software and Instrument Control System Hardware 

The detector is controlled using software provided by Teledyne. Teledyne uses a third party USB 
driver (Bitwise Systems QuickUSB) and their own hardware abstraction layer (HAL) to translate 
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application commands into driver specific commands, hiding the USB driver details. In turn, a 
Teledyne supplied COM DLL objeet uses a .NET web serviee to eommunicate with the HAL. 
This COM DLL object allows third party software to interfaee with the HAL. Teledyne provide an 
example front-end interface to these serviees written in IDL. This interfaee has an inbuilt soeket 
server, enabling eontrol over TCP. For this interfaee to work, a speeial version of the assembly 
eode (HxRG) must be uploaded to the SIDECAR. This eode allows the user to ehange eertain 
exposure options “on-the-fly” without having to reupload the assembly code. Sueh options inelude 
the number of deteetor outputs used during readout, preamp gain, windowing and the exposure 
mode (viz. UTR sampling, fowler sampling) together with their assoeiated parameters. 

As the USB driver is Windows only, 10:1 uses two ICCs: one running Windows and one run¬ 
ning a recent version of CentOS, a stable and predietable derivative of Red Hat Enterprise Einux. 
The windows maehine is used to host the IDE soeket server that eommunieates with the deteetor. 
The linux maehine hosts a remotely accessible instrument agnostic interface whieh eommunieates 
with the soeket server. This interfaee uses a eommand set that is homogeneous to the other in¬ 
struments on the ET,^ allowing for easy integration with the teleseope’s Robotie Control System 
(RCS). The linux maehine also deals with other aspeets of data management, sueh as renaming 
output files in aceordance with ET eonventions and adding the appropriate FITS headers relating 
to the teleseope and observing eonditions. 

2.3.1 A Note Regarding Virtualisation 

It has already been shown by RATIRp^ that this two maehine setup lends itself well to virtualisa¬ 
tion. For 10:1, therefore, these maehines were initially hosted as guests under a linux host running 
VirtualBox, an open souree virtualisation environment. A separate shared storage space within this 
environment was provided to both of the guests in order to store data taken with the instrument. 
However, two key problems were eneountered during testing of this mode of operation: i) sig¬ 
nificant overhead on performing header keyword write/update routines provided by the CFITSIO 
library, and ii) seemingly random FIFO overflow errors during some exposure sequenees. 

The first of these issues arose from the need to append and update headers to the files outputted 
by the IDE software. Quantitatively speaking, it took approximately two seeonds per file to proeess 
the required 72 header cards. Given that the proeess would return almost instantaneously on the 
native ext3 file system, this overhead is most likely a eaehing issue relating to the virtual filesystem 
(VFS) VirtualBox uses for shared folders. The problem was also found, albeit to a lesser extent, 
when using VMWare - a proprietary alternative to VirtualBox. Even for exposure sequenees eon- 
taining a small number of images, this delay beeomes a signifieant eontributor to the instrument 
usage overhead. 

The FIFO overflow errors are more problematic as they eause an abort of the exposure se- 
quenee when they oeeur. This problem is not exelusive to a virtualised setup, but is thought to 
be exaeerbated by the additional overhead ineurred by writing to the VFS. Other solutions were 
tested, ineluding writing to a Network File System (NFS) share located elsewhere on the network, 
but this did not address the underlying eause and the problem persisted with a similar frequency of 
failure. 

To mitigate both of these issues and improve reliability, two separate physieal maehines are 
now used and the files are shared between them both through an NFS share hosted on the windows 
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machineQ. In this manner, the files are written loeally to the drive’s native filesystem, now hosted on 
a solid state disk with mueh greater Input/Output Operations Per Seeond (lOPS) than the previous 
spinning hard drive. The header write/update overhead has notably deereased by roughly an order 
of magnitude and, to this date (Deeember 2015), no FIFO overflow issues have been seen to oeeur. 

2.4 Electrical 

A broad overview of the eleetrieal and optieal eonneetions both inside and outside the eryostat is 
shown in Figure [6l The data lines for the JADE2 are eoupled to a fibre optie transeeiver, allowing 
the transfer of data from the teleseope to an ICC whieh is raek-mounted in a room adjoined to the 
teleseope enelosure. 


TELESCOPE PORT 


PANEL (telescope MOUNT) RACK 



Fig 6 Wiring schematic. REX and LEX refer to the remote and local fibre optic extender units respectively. TLAN 
refers to the telescope ethernet network. 


2.4.1 Noise and Grounding Considerations 

To minimise voltage ripple, two linear power supply units (PSUs) provide power separately to both 
the SIDECAR and the USB seetion of the JADE2 eireuit. A switehed-mode PSU was initially used 
to power the JADE2 USB, but was found to introduee notieable mains piekup noise (see Eigure 
|7]). To avoid ground loops, all internal analog/digital grounds from both the SIDECAR and JADE2 

thttp://www.hanewin.net/nfs-e.htm 
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boards (apart from the JADE2 USB) are tied to a single point ground provided by the SIDECAR 
PSU (see Ref. [8] for speeifies). A replaeement PSU was soureed for the SIDECAR after the noise 
from the first unit was found to eorrelate with increasing temperature, producing vertically banded 
frames upon readout. 
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Fig 7 Noise power spectrum before and after both the faulty SIDECAR PSU was replaced and the switched-mode 
PSU powering the JADE2 USB was changed. The spike at 50Hz in the before PSU replacement plot is mains pickup, 
most likely the result of using a switched-mode PSU. 


3 Detector Characterisation 

The following data was taken in UTR, Correlated Double Sampling (CDS) and Eowler observing 
modes. To minimise analog-to-digital converter (ADC) readout noise, the detector was operated 
at the slower pixel clocking speed of lOOkHz using the Successive Approximation ADCs with a 
resolution of 16 bits. The output signals were measured differentially against a common reference, 
InPCommon, from the detector. The temperature of the array was held at 85K. 

3.1 Dark Current 

During early testing of the array, a blank was used in place of a filter to enable true darks to be 
obtained. However, a combination of reset anomaly and an elevated charge accumulation rate 
made these frames unusable for the purpose of measuring dark current. Without empirical data, 
we must rely on the manufacturer’s specification. Teledyne quote a median dark current for this 
array of less than 0.01e-/s at 120K. The magnitude of this figure is corroborated by RATIR who, 
with the same model array, were unable to measure any dark current for 1000s CDS exposures 
operating at a detector temperature of ~60K.l^ Although we now cannot remove this component 
independently, the dark current, and indeed any other 2D imprint on the data (such as thermal 
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glow from the telescope), is typically removed during the sky subtraction phase of data processing. 
However, if the pipeline-generated sky is a poor representation of the true sky - as will be the 
case if the observation is of a sufficiently extended source - then the observer will be required to 
generate their own offset sky and do this step of the reduction process themselves. 

3.1.1 Reset Anomaly 

In all of the UTR sequences taken with the blank, an appreciable early time nonlinearity was found 
to affect all pixels in the array. This nonlinearity manifested as an additive component to the 
integrated signal, and took an approximately exponential functional form. The problem, known 
as “reset anomaly”, has been observed previously in hybridised HgCdTe detectors and is thought 
to be due to charging effects in the readout integrated circuit (ROIC)/detector.^ Reset anomaly is 
normally only observed to affect a fraction of the pixels in the array, but the effect here presented 
in all pixels.^ Neglecting to take this effect into account leads to systematic overestimation when 
determining dark current. 

A solution has previously been put forward to allow determination of a more accurate value for 
dark current when reset anomaly is present.^ This solution involves determining the coefficients 
of a four parameter fit to the integrated signal as a function of time. These parameters describe 
both exponential and linear components, the latter of which can be used to calculate the dark 
current in an unbiased manner without having to arbitrarily estimate where the contribution of 
the exponential component becomes negligible. However, it is evident from the magnitude of the 
charge accumulation rates in Figure [5] that the dark frames obtained for this analysis contain an 
additional diffuse component that uniformly illuminates all pixels in the array. This component is 
thought to be due to radiation originating from the blank and the surrounding filter housing. As 
the contribution from this component can be subtracted off during the data reduction process, this 
problem has not been investigated further and consequently no direct empirical measurement of 
dark current is possible. 

Additionally, the form and magnitude of the reset anomaly for this detector has been observed 
to change with time. Although the effect was not observable in any of the test data used in the 
following sections, it has presented again in UTR sequences taken at several dates following a 
recent thermal cycle. Plots of the mean count as a function of time are shown in Figure [U The 
detector had settled at its operating temperature. 

To ascertain if this effect was primarily additive (and thus fully correctable during data reduc¬ 
tion), the photometric linearity tests discussed in ^5.II were performed on datasets both with and 
without reset anomaly. No measurable error on the former dataset could be detected. 

3.2 Gain and Read Noise 

In the absence of true dark frames, the Photon Transfer MethocP^ was used to measure both read 
noise and gain for a range of preamp gains from 12dB to 21dB in 3dB steps. For each preamp gain, 
the reference voltage Vrefmain (SIDECAR register 0x602c) was adjusted so that the reference 
pixel signal from the detector was in range of the ADCs. The array was partitioned into 50x50 
pixel windows and the gain/read noise calculated for each. A histogram of the results for each gain 
setting can be seen in Figure l9l and [TOl 

For the four preamp gains of 12, 15, 18 and 21 dB, the four conversion gains attainable are 
3.05, 2.12, 1.5 and 1.06 e-/ADU respectively. 
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Fig 8 An illustration of how the reset anomaly is changing with time. Acquisition dates for each of the sequences are 
shown in the legend. In the most recently taken sequence (2015-06-09), the nonlinear jump in magnitude during the 
early part of the ramp is seen to be considerably smaller. The later part of sequences 2015-06-03 and 2015-06-09 are 
also more linear than 2015-05-26 after ~5s. 

The reference subtracted read noise calculated in Figure [TO] is the CDS read noise, which is 
greater than read noise for a single frame by a factor of \/2. By observing in Fowler mode, the 
read noise can be further decreased by increasing the number of pairs in the observing sequence; 
this is illustrated in ^3.5[ For the four preamp gains of 12, 15, 18 and 21 dB, the four values of 
CDS read noise attainable are 32.4, 26.5, 21.8 and 19 e- respectively. 

3.3 Full Well Depth and Nonlinearity 

Full well depth (FWD) was measured by taking long UTR sequences that either saturated the 
detector or the ADCs. The results for each gain setting can be seen in Figure [H] To match 
Teledyne’s definition, FWD is hereinafter defined as the count at which nonlinearity >5%. 

For the four preamp gains of 12, 15, 18 and 21 dB, the four FWDs attainable without correction 
are 102,72, 69 and 57 ke- respectively. Given that a higher preamp gain both reduces the read noise 
and full well depth, a compromise has had to be sought by distillation of the preceding analysis to 
the following statements: 

- The 12dB setting underutilises the full range of the ADC and has the highest read noise. 

- The 21dB setting has a conversion gain of ~unity, and is thus only capable of a FWD of 
< 2^® e-. Given that a typical reduction sequence for the data is subject to both the removal 
of a DC bias level and CDS, this figure is actually an overestimate. Empirically this value 
has been found to be ~57ke-. 
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Fig 9 Conversion gain for the different preamp settings. Conversion gain decreases as preamp gain increases. 


- Both 15dB and 18dB settings are viable. The 15dB setting has the advantage of only ~3ke- 
larger FWD, but 18dB has ~4.7e- less read noise. 

Favouring the reduetion in read noise over inereasing the FWD, the preamp gain has been set 
to 18dB. 


3.4 Nonlinearity Correction 

To aeeount for nonlinearity, a series of ealibration files were eonstrueted using UTR sequenees. 
This was done for only the 18dB setting. For eaeh pixel, the integration time was plotted against 
the eounts observed. A linear fit was then made to the early part of this sequenee when the eharge 
aeeumulation is assumed to be linear (no reset anomaly was present), and a third order polynomial 
fit made between the values extrapolated from the linear fit - the “aetual” eounts - and the observed 
eounts. Per-pixel eorreetion eoeffieients were then ealeulated. 

A eorreetion derived from these eoeffieients has been applied to a separately obtained UTR 
sequenee and the residual nonlinearity ealeulated. The result is shown in Figure [T^ 
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Fig 10 CDS read noise for the different preamp settings. The read noise was converted into units of electrons by using 
the conversion gains determined previously. Read noise decreases as preamp gain increases. 


3.5 Number of Fowler Pairs 

For Fowler sampling, read noise is dependent on the number of pairs taken. The greater the number 
of pairs, the smaller the read noise that is attainable through averaging. This is shown in Figure [T3l 
This redueed read noise comes at the cost of increased observational overhead. For each additional 
pair, an extra read at the end of the sequence is required. Operating in lOOkHz 32 output mode, 
a single extra read corresponds to an additional ~1.3s. Increasing the number of pairs also has 
implications for data transfer from site and the resulting processing overhead. 

3.6 Summary 

A summary of the preceding analysis is shown in Table [B 

4 Data Pipelining 

When used through the telescope’s normal robotic interface, the instrument automatically se¬ 
quences for the observer a dithered series of exposures to allow sky subtraction. 
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Fig 11 Measurements of full well depth. The results for the different preamp settings are shown on different rows. 
The left panes show how the mean count increases for a region of the array that was free from cosmetic defects. The 
middle panes illustrate how the CDS mean count deviates from a linear fit to early integration times. The CDS mean 
count is taken as the difference between the first and last reads. The right panes show the percentage nonlinearity as a 
function of the mean count. 


Table 1 A summary of lOT’s characterisation tests. 


Preamp Gain 

12 

15 

18 

21 

Conversion Gain (e-/ADU) 

3.05 

2.12 

1.50 

1.06 

CDS Read Noise (e-) 

32.4 

26.5 

21.8 

19 

Uncorrected FWD 5 % (ke-) 

102 

72 

69 

57 

Corrected FWD 5% (ke-) 

- 

- 

93 

- 

Vrefmain (V) 

1.71 

1.26 

1.16 

1.09 


An autonomous pipeline has been developed to reduee data taken using this standardised ob¬ 
serving sequenee. This pipeline performs referenee subtraetion, frame eombination, non-linearity 
eorreetion, flatfielding, bad pixel masking, sky subtraetion, registration and staeking. It is written 
in Python using PyFits and the third-party modules aiipy@ and statsmodeisH. An example re- 

^http://obswww.unige.ch/^tewes/alipy/ 

§https://github.com/statsmodels/statsmodels/ 
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Fig 12 Nonlinearity before and after correction for an UTR sequence at a preamp gain of 18dB. 

duction sequence is shown graphically in Figure [HI A brief description of the reduction recipe 
follows. 

4.1 Reference Pixel Subtraction 

The active pixels of the H2RG array are bordered by a total of 8 columns and 8 rows of inactive 
pixels. Although these inactive pixels are not bonded to a photodiode, they are clocked during 
readout giving a means by which to reduce the common-mode noise present in active pixels. The 
simplest implementation of reference subtraction is to take an average of all of the inactive pixels 
and subtract this value from that of all active pixels. However, doing so does not take into account 
the cross-array, column-to-column and output-to-output drifts. 

In an attempt to correct for these nuances, several different reference pixel subtraction schemes 
were tested. These schemes included, in varying combinations, single and per-output subtrac¬ 
tion, fitting a linear slope to the upper and lower reference pixel regions (ramp) and distinguishing 
between odd-and-even columns (odd-and-even). For those schemes employing per-output subtrac¬ 
tion, the array was split into 32 sections (i.e. the full 2048 pixel width of the array is split into 
sections of width 64 pixels, and each section is read through a different output amplifier). For 
those schemes employing ramp subtraction, the offset to be subtracted is a function of vertical po¬ 
sition of the pixel in the array and was calculated on a per-row basis. A scheme to remove residual 
vertical banding by smoothing the columnar reference pixels is also currently being investigated. 
A comparison of the resulting read noise for each of the different methods is shown in Tabled 
In practice there is no measurable difference between subtracting a single per-output constant 
and fitting a ramp to the data. Distinguishing between odd and even columns also has no impact 
on the read noise. The pipeline currently uses per-output ramp -i- odd-even reference subtraction. 
More complex subtraction schemes do exist that make better use of these reference pixels, but 
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Fig 13 Average read noise as a function of number of fowler pairs for a preamp gain of 18dB. 

Table 2 Read noise for a variety of reference pixel subtraction methods. These values were calculated with a preamp 
gain of 18dB. 


Method 

Read Noise (e-) 

None 

24.4 

Single frame constant 

25.2 

Per-output constant 

21.8 

Per-output constant + odd-and-even 

21.8 

Per-output ramp 

21.8 

Per-output ramp -i- odd-and-even 

21.8 

Per-output ramp -i- odd-and-even -i- vertical banding removal 

21.1 


these sehemes have not been explored as they require alteration of the eloeking pattern for readout 
of both referenee and aetive pixels.^ 

4.1.1 Correlated Noise 

Due to temporal drifts in the bias voltages supplied to the deteetor by the SIDECAR, there is the 
possibility of correlated noise between amplifiers. To quantify the extent this noise, we separated 
a short exposure image into 32 strips of data corresponding to one strip per amplifier. Each strip 
was then differenced with each other strip, and the standard deviation of this difference compared 
with the theoretical value expected for uncorrelated noise, given by taking the standard deviations 
of the two strips and adding them in quadrature. The distribution of expected theoretical noise 
relative to that which was measured is shown in Eigure [JS] There is evidence of a comparatively 
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Fig 14 Example intermediate and end products from a typical reduction sequence using two pairs and N dithers. To 
make the stack, the dithers are averaged and scaled up by a factor of N. 

low level (<3 ADU) of correlated noise (cf. ~24 ADU from Table O. This issue may be further 
investigated and addressed in a later version of the pipeline. 

4.2 Frame Combination 

Both CDS and Fowler sampling are supported by the pipeline. As discussed in 1 13.51 Fowler sam¬ 
pling with two pairs is used for routine robotic operations. 

4.3 Nonlinearity Correction 

The calibration files discussed in 1 13.41 are used to correct for nonlinearity. Each order of the correc¬ 
tion is stored in a separate extension of a single multi-extension file, allowing for easy application 
by array multiplication. 

A caveat exists to performing nonlinearity correction when the combination scheme is either 
CDS or Fowler sampling. In order for the nonlinearity correction factor to be calculated using an 
accurate reflection of the true charge stored in a pixel, it is required to first estimate the flux that 
is “lost” from the first readout due to taking the difference between pairs of frames. The timescale 


17 









































60 


50 


40 


E 30 


20 - 


10 


ill 


-4.0 -3.5 


3.0 -2.5 -2.0 -1.5 -1.0 -0.5 

Deviation from uncorreiated noise (ADU) 


0.0 


0.5 


Fig 15 A histogram of the deviations between the expected theoretical noise of two differenced strips and the mea¬ 
sured value. A negative value indicates that the measured noise was less than theoretically expected, indicating that 
correlated noise may be present. The bimodal form of the distribution indicates that only certain pairs of amplifiers 
are affected. 

over which this flux is lost is equivalent to a multiple of the readout time. Calculating the rate 
of flux accumulation between the first and second frames of the first pair and factoring this into 
the linearity correction ensures that the correction is applied to the true value of pixel charge. If 
this lost flux is not taken into account, brighter targets (those with a higher rate of flux) will be 
systematically undercorrected. 

4.4 Flatfielding 

A variety of flatfielding techniques were considered but were deemed too problematic to obtain 
without a cooled blank. Without a blank, it is not easy to disentangle the QE of the detector and 
the dark current component. Doing so requires exposures of the same length but with different 
background brightnesses. Given this limitation, we have decided to use twilight difference images 
to construct the Hatfield. There will inevitably be some discrepancy introduced from using twilight 
frames with a much higher colour temperature than typical of the sources being observed, but given 
that other schemes are not accessible to us presently, this is an unavoidable caveat. It should be 
noted that other flatfielding schemes, such as those using stacks of the night sky, may be ideal for 
flatfielding the night sky but would not necessarily match the colour of observed targets either. 

To construct our master Hatfield, a library of difference images was assembled using tracked 
and dithered twilight images of a blank field. Each image was processed through the normal 
pipeline as far as and including the Eowler combination. Difference images were then created 
using pairs of bright and faint twilight images of the same integration time. This removes all 
additive noise components such as overall dark current (expected to be negligible), dark current 
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from a few hot pixels, reset anomaly and thermal baekground glow from the teleseope itself - 
typically of order IkADU - assuming that this glow is constant on short timescales (any change in 
this background would imprint on the flat). These difference images are level matched, combined 
with deviant pixel rejection and normalised to create the master Hatfield. 

4.5 Bad Pixel Masking 

By applying constraints to both flatfielding coefficients and nonlinearity residuals, pixel maps have 
been constructed to filter out bad pixels which either have a QE of <0.35 of the average and/or 
median residual nonlinearity of >0.8%. The distribution of these two parameters is shown in 
Figure [m 




Fig 16 Distribution of pixels as a function of flatfielding coefficient value (top) and residual nonlinearity (bottom). 
The dashed lines represent the lower (flatfielding) and upper (residual nonlinearity) limits set for a pixel to be flagged 
as “bad”. 

4.6 Sky Subtraction 

Sky subtraction is currently achieved by construction of a median stack using all dither positions 
except the one for which the stack is being generated (“peers-only”). At the cost of increased sky 
noise, peers-only combination reduces the likelihood of systematically overestimating the central 
bright regions of each stellar PSF when using lesser robust estimators such as the median. An 
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iterative sigma elip is then applied to determine and mateh the baekground levels of both the dither 
position being eonsidered and the staek. The baekground-level-matehed staek is then subtraeted 
off. 

The sky subtraetion stage of the pipeline is something of a misnomer, insomueh that this stage 
also removes the eontribution from both dark eurrent and thermal glow. 

4.6.1 Unbiased Estimation of the Sky Value 

Other more robust M-Estimators were eonsidered and tested, ineluding Tukey’s Biweight and Hu¬ 
ber’s T. Both of these algorithms use an iterative maximum likelihood solver that strives to seleet 
the “average” value that is most eonsistent with the sample values being drawn from a symmetri- 
eal, normal distribution. In the eontext of this partieular applieation, finding the average value is 
equivalent to estimating the value of the sky at a given pixel. Like the mean, all data are eonsidered 
in the final average but instead of eonsidering eaeh data point with equal weight, the data are mul¬ 
tiplied by a weighting funetion that down-ranks points that lie further from the determined eentre 
of the distribution. 

Initial testing of Tukey’s Biweight looked to be the most promising alternative, proving mueh 
more robust to outliers than using a median or Huber’s T. However, when applied to real data it 
was found that while the algorithm provides a better esimator of the sky baekground at the eentre 
of the star PSFs the outer haloes of the stars were oversubtraeted resulting in ringed artefaets (see 
Figure [TT]). 



Fig 17 Illustration of residual artifacts post sky subtraction formed when using both a simple median and Tukey’s 
Biweight to generate the stacked sky. The example is from a 4 dither pattern. Panels a) and b) are stacks generated 
using a median. Panels c) and d) are stacks generated using Tukey’s Biweight. Panel a) and c) were generated using 
all frames. Panels b) and d) were generated using the peers-only scheme described in ^4.61 The sky noise in panels 
a), b), c) and d) is 37.7, 51.1, 36.6 and 49.4 ADU respectively. 
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With careful consideration of aperture size, Tukey’s Biweight yields the best absolute photom¬ 
etry. However, the artefacts created were thought to unnecessarily introduce confusion for users 
who were unfamiliar with this algorithm. This is especially true when one considers that the differ¬ 
ence between the various subtraction schemes is significantly less noticable when an offset pattern 
utilising a greater number of dithers is selected. Consequently, the use of Tukey’s Biweight was 
rejected for common use in favour of a simple median. 

4.7 Image Registration and Stacking 

The pipeline uses the aiipy module to find the per-frame tranformation matrix with components 
of translation, rotation and scale. The first dither position is selected as a reference by which all 
the other dither positions are registered to. 

5 Photometric Tests of the Pipeline and Instrument Performance 

The best test of the reduction pipeline and real world instrument accuracy is the consistency and 
repeatability of photometric measurements on sky. For this analysis, four fields were selected and 
observed with various integration times ranging from 6 to 60 seconds. Though the selected fields 
were centred on stars from the ’UKIRT Fundamental and Extended Lists’]^ all 7 < H < 16 stars 
in the 2MASS Point Source Catalog were included in the analysis unless otherwise stated. 

5.1 Linearity and Sky Brightness 

To assess photometric linearity, observations of these fields were run through the 10:1 pipeline and 
sources extracted using the Source Extractor (SExtractoi^ software with several different circular 
apertures between six and twenty pixels diameter. Counts were converted to photo-electrons per 
second using a gain of 1.5e/ADU and instrumental magnitudes calculated. EigurefTSl shows a col¬ 
lation of these instrumental magnitudes on three different nights for the field surrounding standard 
star ES150, plotted against 2MASS PSC H-band magnitudes. The gradient of the fit to the data is 
consistent with unity, signifying the detector has responded linearly in these observations. 

Since the detector red cut-off lies within the H-band it is to be expected that there will be a 
colour transformation between lOLI H-band and other instruments. In practise, the transformation 
is not large for typical stellar sources though it could become significant for any particular source 
with strong spectral features between 1700 and ISOOnm. Eigure [TS] also shows the colour trans¬ 
formation with respect to 2MASS colour has a slope of -0.093 x (J-H) for all stars in the ES150 
field. 

The magnitude of the sky background was also calculated in the same analysis. We report 
a value of 13.34±0.08 mags/arcsec^ - about a half a magnitude brighter than the values quoted 
by both TNG0 and INT/IKlfi]. This is not surprising given that no NIR optimisation has been 
made to either the telescope’s infrastructure or environment. Indeed, there is a difference of ~0.3 
mags/arcsec^ between the centre of the field and the outside, presumably the result of residual 
thermal glow. 

^http;//www. tng.iac.es/instruments/nics/imaging.html 

II http;//www.ing.iac.es/Astronomy/observing/conditions/skybr/skybr.html 
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Fig 18 Top; Instrumental magnitude of sources calculated using frames taken with 10:1 compared with corresponding 
values taken from the 2MASS catalogue. The overplotted straight line is not a fit, but simply a gradient of unity which 
would result from using a perfectly linear detector. Free fits to the data from different fields yielded gradients in the 
range 0.996 to 1.004. Bottom: Colour transformation with respect to 2MASS colour. 

5.2 Consistency ofinvidual Frames and Stacked Products 

To check that the aligned, co-added data are stacked correctly in the pipeline, photometry was 
extracted independently for the same stars from each individual frame and from the final stacked 
data product. The zeropoints derived for all stars in two image fields are plotted against exposure 
time in Figure [19] for integrations ranging from 6 to 60 seconds and their associated five-frame- 
stacks. Photometry was taken from the sky-subtracted pipeline images with faint and saturated 
stars excluded. 

5.3 Sky Subtraction 

To test the pipeline sky subtraction and ensure that star residuals were not imprinted on the derived 
sky image, photometry was extracted from both the sky-subtracted image (FITS extension IM-SS) 
and the non-sky-subtracted data product (FITS extension IM-NONSS). Photometry was performed 
in SExtractor using both simple circular apertures and the ’’AUTO Kron-like elliptical aperture” 
mode. The two extraction modes gave consistent results. For the sky-subtracted image no sky 
subtraction was performed in SExtractor. Eor the non-sky-subtracted image, sky subtraction was 
performed using the package’s default filtered sky background algorithms. Eigure|2Q|compares the 
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Fig 19 Zeropoints calculated from both individual exposures and the stack as a function of exposure time. The 
zeropoint is defined to be 1 ADU/sec referenced to the 2MASS PSC natural system and with no colour transformations. 

ensemble results for six different integration times on two different star fields. The derived ZP with 
respect to 2MASS is plotted both against total counts in the extraction aperture and against peak 
brightness at the centre of the star. Results are seen to be consistent for the two datasets and the 
scatter significantly reduced by using the pipeline’s automated sky subtraction. This is expected 
since the pipeline makes use of multiple stacked images (see ^4.61) to derive a clean sky image with 
all the stars and other noise sources such as hot pixels removed. 

5.4 Persistence 

Image persistence from bright sources into subsequent images is detectable. A dithered sequence 
of images centred on a heavily saturated star with estimated peak counts of 1.3 x lO^ADU shows 
a persistent image of peak 350ADU one minute after the initial exposure, fading to 170ADU 
two minutes after the exposure but undetectable after four minutes. Persistence is therefore not 
normally expected to be a significant problem and even one minute immediately after the saturation 
event the residual is less than 0.03 per cent, likely smaller than many other sources of photometric 
error. It is noted that using this analysis of persistence doesn’t entirely probe the mechanisms 
whereby charge is trapped, which has been shown to be dependent more strongly on exposure 
cadence rather than intensity.!^ However these tests were performed with 60 second integrations, 
the longest recommended given the sky brightness and thermal background, thus representing an 
upper limit of persistence expected in normal science operations. 

6 Conclusions 

In this paper, the mechanical, electronic and cryogenic aspects of lOT’s development have been 
presented. This was followed by a discussion of results derived from characterisation tests, includ- 
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Fig 20 Peak count rates and total counts as a function of zeropoint calculated with both default pipeline sky subtraction 
and without (instead using SExtractor sky subtraction). The striations in the lower panel are the result of using several 
exposures of the same field. Each group of points represents many observations of the same star. These groups form 
tilted lines because the 2MASS comparison reference which the observed counts are compared against is a single 
catalogue value whereas there is numerical scatter between our multiple observations. 

ing measurements of read noise, eonversion gain, full well depth and linearity. These measure¬ 
ments predieated the seleetion of 18dB as the most appropriate preamp gain, with eorresponding 
parameters given in the 18dB eolumn of Table [B An overview of the data pipelining proeess was 
given and the results of photometrie tests eondueted on on-sky, pipelined proeessed data were pre¬ 
sented. These tests have proved that the pipeline produees eonsistent and repeatable data produets 
in terms of photometrie linearity, both in sky-subtraeted and non-sky-subtraeted frames. 

Observers interested in solieiting a proposal to use 10:1 are eneouraged to apply through the 
website at http://teleseope.astro.ljmu.ae.uk/PropInst/Call/, To assess programme suitability, an ex¬ 
posure time ealeulator is available at http://teleseope.astro.ljmu.ae.uk/TelInst/eale/. 
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